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Abstract 

Organic Zener diodes with a precisely adjustable reverse breakdown from -3 V to -15 V 
without any influence on the forward current- voltage curve are realized. This is accomplished 
by controlling the width of the charge depletion zone in a pin-diode with an accuracy of one 
nanometer independently of the doping concentration and the thickness of the intrinsic layer. 
The breakdown effect with its exponential current voltage behavior and a weak temperature 
dependence is explained by a tunneling mechanism across the HOMO-LUMO gap of neigh- 
boring molecules. The experimental data are confirmed by a minimal Hamiltonian model 
approach, including coherent tunneling and incoherent hopping processes as possible charge 
transport pathways through the effective device region. 
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In the last few years, the physical properties of organic semiconductors have been studied 
with considerable detail. Despite the fact that there is still a limited understanding of the relevant 
microscopic charge transport scenarios in such systems, various organic electronic devices such 
as efficient organic light-emitting diodes (OLEDs) 1 ' 2 , organic memory cells 3 , organic field-effect 
transistors 4 , vertical Modes 5 , and high-frequency diodes 6 have been successfully demonstrated. 

The main difference between inorganic and organic semiconductors is that the latter usually do 
not display (at room temperature) band-like transport through delocalized states, but rather exhibit 
hopping of carriers between spatially localized states. This can either happen as electron transport 
through the lowest unoccupied molecular orbitals (LUMO) or hole transport through the highest 
occupied orbitals (HOMO). Typically, this results in mobilities three or more orders of magnitude 
below those of band transport and an entirely different temperature dependence, i.e. an increase of 
mobility with increasing temperature can be obtained. 

In inorganic semiconductors, there are a number of basic transport effects which involve more 
than one band, e.g. interband (Zener) tunneling or avalanche breakthrough. 7 ' 8 While these effects 
are now understood in detail for inorganic semiconductors, they have not been experimentally 
studied in detail for organic semiconductors so far. Due to the structural complexity of organic 
semiconductors, the theoretical study of their charge transport properties is also quite challenging, 
so that theoretical studies are usually limited to structurally ordered systems. 9-12 A key electronic 
device based on the Zener effect, the Zener diode, is crucial for basic electrical circuit requirements, 
such as voltage- and temperature stabilization and over- voltage protection. An organic Zener diode 
with an adjustable breakdown is still missing. Likewise Zener diodes are key elements for passive 
matrix memory arrays, 13 where Zener diodes are essential to prevent parasitic current flow through 
non-selected crosspoints. The required parameters for the Zener diode are given by the read-, write- 
and erase voltages of the memory. 14 

In the Letter, we investigate the HOMO-LUMO Zener tunneling in an organic semiconductor 
system by detailed measurements of the reverse currents in organic pin-diodes. In particular, we 
provide evidence for the tunneling nature of the reverse currents by showing a strong influence 
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of the effective tunnel barrier and a very weak temperature dependence. To get a first insight 
into the possible charge transport mechanisms, we complement the experimental investigations 
by an effective low-dimensional model which is then treated in two limiting physical situations: 
incoherent hopping transport and coherent tunneling. Our results suggest that both mechanisms 
may exist in the system under study in dependence of the sign of the applied voltage. 

Experimental part. We address Zener tunneling in the reverse current regime of organic pin 
diodes. The breakthrough properties of such devices have only been rarely investigated, see e.g. 15 
The advantage of the pm-diode concept 16 is that the geometry and the potentials within the semi- 
conductor layers can be controlled very precisely, thus allowing to study the tunneling effects, 
which sensitively depend on minute variations of these parameters, in a very controlled man- 
ner. Thus, the tunneling current is independently tunable by two parameters: the thickness of 
the intrinsic interlay er (IIL), and the doping concentration of the electron and hole transport layers 
(ETL/HTL). Doping changes the charge carrier concentration 17 and width of depletion regions. 18 
Furthermore, one obtains a barrier-free transport until the charge carriers reach the IIL. 19 In con- 
trast, the IIL acts as an additional depletion layer and the voltage drop over the pin -junction mainly 
occurs in this intrinsic layer. 19 A comparative study of organic pn- and p?>z -junctions comprising 
doped organic layers would be desirable. However, as we will show in this letter the typical width 
of the native depletion zone in the doped layers is less than 3 nm. Therefore, the current in a in- 
junction would be caused not only by drift and diffusion but also by a tunneling mechanism. In 
this way the backward current should not be tunable independently of the forward direction. 

The structures investigated (see the inset of [figure] [1][] lb) consist of three organic layers 
sandwiched between two aluminium contacts, which define the active area of 6.38 mm 2 . The inner 
structure of the organic layers follows the />m-concept with a highly p-type doped HTL, a highly 
n-type doped ETL, and an intrinsic layer in between (see Supporting Information for details), [fig- 
ure][l][]l displays the current- voltage (I-V) characteristics of such a device. In forward direction 
the expected diode-like behavior 16 is observed. The interlayer thickness does not affect the onset 
and the slope of the current rise in forward direction, which appears at 2.3 V. In contrast, there is a 
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Figure 1: Current voltage characteristics for different interlayer thicknesses (b) Capacitance volt- 
age characteristics (measured at 1 kHz) of the samples shown in [figure] [1][] la. The inset is 
a sketch of the device structure and the layers are listed as they are evaporated on a cleaned 
glass substrate (1): Al(100nm)(2) /MeO-TPD:F4TCNQ[4wt%](50nm)(3) /Balq2:NPB[50wt%](4) 
/Bphen:Cs(50nm)(5) /Al(100nm)(6) 
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significant influence on the reverse current, which can be seen in the loss of the rectifying property 
at a certain voltage depending on the interlayer thickness. The onset of the reverse breakdown 
shifts in a controlled manner to higher negative voltages (from -3 V to -15 V) by increasing the 
interlayer thickness (see [figure] [1][] la): we observe a shift of the reverse breakdown voltage of 
approximately 1 V by each additional nanometer of IIL thickness. Current densities of more than 
300 mA/cm 2 in reverse direction can be achieved without limiting reversibility or reproducibility. 
Obviously, the large reverse voltage region is dominated by a purely exponential rise of the cur- 
rent, which provides a strong hint for a tunneling like mechanism. Also a flat part of the I-V curve 
can be observed in reverse direction. Since this is more pronounced for thicker interlayers it is 
considered to be related to a leakage currents. Hence, a larger field and a larger voltage are needed 
to be in the tunneling regime for thicker interlayers. The limitation of the transport layers and the 
injection gives rise to the flattening of the I-V curve for larger negative voltages. 

To further understand these observations, we perform capacitance-voltage (C-V) measure- 
ments, allowing to study the voltage-dependent effects of charge depletion and accumulation in 
the IIL and allow correlating layer thicknesses to capacitance values. In our devices, the transport 
layers have a larger thickness and a better conductivity in comparison to the IIL which is thus 
the dominating 7?C-element in the stack. Capacitance frequency (C-f) measurements are done in 
order to prove independently the capacitance calculated by the assumption of one dominant RC- 
circuit (see Supporting Information). Thus, it is checked that the calculated capacitance at V 
corresponds to the assumed interlayer thickness, which means that the IIL is still depleted even 
for such thin films and no diffusion of dopants occurs. The C-V curves shown in [figure] [1][] lb 
obey the Mott-Schottky rule as long as no significant current appears (between V and -4 V, see 
Supporting Information). In the forward voltage direction thus, a sudden capacitance drop appears 
independently of the IIL thickness at 2.3 V, where the geometrical capacitance dominates the C-V 
curves since the flat-band condition is reached. The small drop of capacitance for a small forward 
bias is related to the ETL system of cesium doped BPhen. The significant number of trap states in 
this material which will be filled under forward voltage conditions causes the drop in capacitance. 
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As one can see e.g. in [figure] [3] [] 3b this effect does not appear in doped RE68. However, in re- 
verse direction, a weak decrease of the capacitance is observed as predicted by the Mott-Schottky 
rule. This corresponds to a charge depletion at the HTL/IIL and ETL/IIL interfaces. Besides this, 
a drop of the capacitance is obtained as expected from the I-V curves for certain reverse voltages 
depending on the IIL thickness caused by a significant injection into the ILL. 

To demonstrate the control of the reverse behavior and to discuss the possible relevant charge 
transport mechanisms, we plot in [figure] [2] [] 2 the reverse current for different temperatures and 
IIL thicknesses. For decreasing ILL thickness (from 11 nm to 5 nm), a purely exponential current 
dependence is observed (for reverse voltages larger than -3 V), with an increase of the reverse 
current by nearly six orders of magnitude. 

A further strong evidence of tunneling is the observed temperature dependence: In addition 
to the exponential behavior, we find a rather weak influence of the temperature on the reverse 
current. In contrast to a typical activation energy of the forward current (is£? f r =(170±20) meV) 
for such a pm-device, 20 a surprisingly small activation energy of E™ v t =(30±l0) meV is measured 
in the reverse direction, without any voltage dependence (see inset of [figure] [2] [] 2b). This low 
reverse activation energy can be interpreted as being related to HOMO-LUMO tunneling processes 
which become dominant in this voltage region. The difference between this proposed mechanism 
and the effect described by Drechsel et al. 21 can clearly be seen in the weak interlayer thickness 
dependence of the backward current reported by Drechsel. This is mainly due to a change in 
the injection conditions for majority charge carriers at the schottky contact and not related to a 
tunneling of minority charge carriers as reported here. Concerning the activation energy in reverse 
direction it should be mentioned that since we can not neglect the influence of the transport layers 
to the temperature dependence the real value of activation energy remains elusive. However, even if 
the found temperature and interlayer dependences are a strong evidence for a tunneling mechanism 
a second effect, the Avalanche breakdown, could also play a role as a possible transport mechanism. 
Since we obtain an exponential current-voltage characteristic for different interlayer thicknesses 
and for a large voltage region just one kind of these two processes is supposed to be dominant. 
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Figure 2: Current density at a reverse voltage of -5 V for different interlayer thicknesses. Device 
structure as in [figure] [1][] 1 (b) Temperature dependence of the current density in reverse direction 
measured at -3 V and in the exponential rise of the forward direction at 2.6 V. The inset shows 
the temperature dependence for different backward voltages. The devices consist of Al(100 nm)/ 
MeO-TPD:NDP2[4 wt%](50nm)/ Balq2:NPB[50 wt%](7nm)/ Bphen:Cs(50 nm)/ Al(100 nm) on 
a cleaned glass substrate. 
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Figure 3: (a) Current density for different doping concentrations in the HTL. The devices consist 
of Al(lOOnm)/ RE68:NDP2[x wt%](50nm)/ RE68(7nm)/ RE68:NDNl[16wt%](50nm)/ Al(100 
nm) on a cleaned glass substrate. Inset: comparison of different dopand materials regarding 
the I-V performance. The structure consists of a doped HTL(50nm), an intrinsic interlayer 
(Balq2:NPB[50wt%] (7nm)) and a doped ETL(50nm). The colors denote (HTL/ETL): Meo- 
TPD:F4TCNQ[4wt%]/ BPhen:Cs (light green), Meo-TPD:NDP2[4wt%]/ BPhen:Cs (magenta), 
Meo-TPD:NDP2[4wt%]/ BPhen:NDNl[16wt%] (dark blue) (b) Capacitance voltage characteris- 
tics of the samples shown in [figure] [3] []3a 8 



Avalanche and Zener breakdown are distinguishable by the field, which needs to be applied to be 
in the breakdown region. 7 ' 8 Thus, for a certain energy gap of the materials used for the diode Zener 
tunneling should appear for lower fields than Avalanche multiplication. Since we do not obtain a 
second process in the current- voltage curves for higher voltages the Zener effect is supposed to be 
the dominant mechanism. 

If one assumes a transport mechanism similar to inorganic Zener diodes, one can further ex- 
pect a control of the reverse breakdown by the doping concentration of the transport layers, [fig- 
ure] [3] [] 3 shows data a for a diode of RE68 designed as a homojunction: the breakdown voltage 
can be shifted to higher values for lower doping concentrations as expected. The reason is a change 
of the effective thickness of the intrinsic layer by an additional depletion region at the HTL/IIL 
and the IIL/ETL interface. Although the reverse current has a reduced sensitivity to the doping 
concentration than to the IIL thickness, this allows a second independent approach to control the 
breakdown voltage. Furthermore, we also obtain no influence on the forward direction of the diode, 
which means that neither the series resistance of the transport layers nor the flat-band condition 
have been significantly changed. 

C-V spectroscopy investigations (see [figure] [3] [] 3b) once more reveal a Mott-Schottky like 
behavior for different doping concentrations. Obviously, there is an increase of the capacitance for 
a decreasing concentration of dopant molecules. Hence, if the thickness of the depletion layer is 
calculated from the capacitance, one obtains e.g. 8.3 nm for the 4 wt% of dopants and 1 1.3 nm for 
the 0.5 wt% of dopants (including the IIL thickness of 7 nm). This means that the depletion region 
thickness is the sum of the IIL and the interface depletion thickness, being sensitive to the doping 
concentration. We obtain a similar rule as mentioned before in case of the thickness controlled 
reverse current: the reverse breakdown can be shifted by 1 V or the reverse current can be increased 
by one order of magnitude by adding one additional nanometer to the thickness of the depletion 
layer (IIL plus interface depletion). This doping controlled reverse behavior can also be obtained 
for various concentrations of n-type dopants in ETL without any noticeable difference. Instead of 
NDN1 one can use cesium as a freely available n-type dopant (see Supporting Information). 
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Theoretical Modeling. The IIL part of the device is not expected to have an ordered crystalline 
structure, but rather display a high degree of static and dynamic disorder. Since the spatial arrange- 
ment of the molecular building blocks in the IIL is not known a priori, a first-principle description 
of the IIL electronic structure as well as of charge migration 9 " 11 in such a system represents a very 
strong challenge. Hence, to gain a first qualitative insight into the possible charge transport path- 
ways, we have adopted a minimal model approach. Basically, the model consists of an electronic 
ladder with N blocks describing the IIL, each block containing two energy levels ej(V) whose en- 
ergetic position is determined not only by the intrinsic electronic structure of a molecular building 
block but also by the built in field and by the applied voltage (see inset of [figure] [4] [] 4). The 
Hamilton operator for this tight-binding model can be written as: 

H = E %$<y)d + d sj + £ t s ^(d+d sj+1 +h.c.) 

s=H,L j s=H,L j 

+ t%- L Y,(d£ J d Lj+1 + hx.)+tjl- L Y,(d£ j d Lj+ 2 + h.c.) 
j j 

+ £ £ka4 a Cka+ £ £ t kia j :S (c+ a d sj + h.c). (1) 
k,a=l,r k,a=l,r j,s=H,L 

The levels are denoted by s = H (highest occupied molecular orbital, HOMO) and L (lowest- 
unoccupied molecular orbital, LUMO) and are thus representative of the frontier orbitals of a 
given molecular unit. Their energetic position is assumed to be given by the linear relation: 
£j(V) = £ s — Kj + jjeV, where £ s is the bare onsite energy of the molecular orbital (assumed 
to be site- independent for simplicity). The local H-L gap A# -l = £ w — £ L ~ 3.0eV is taken from 
the experimental energy diagram profiles. The built-in electric field is modeled as a linear ramp 
with slope — K, which is supported by the current experiments as well as by photo-electron spec- 
troscopy measurements in organic pm-diodes. 19 The second, third, and fourth summands in Eq. 1 
describe different electronic couplings within the ladder: intra- strand nearest-neighbor couplings 
(second term) as well as nearest- and next-nearest neighbor inter-strand interactions (third and 
fourth terms). If not stated otherwise, we will assume for the sake of simplicity in the calculations 
t H = t L = to and t^~ L and tjf ~ L = rjto, with < r\ < 1. All electronic parameters can in princi- 
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pie be obtained from microscopic first-principle calculations; this goes however beyond the scope 
of the present minimal model approach. The last row in Eq. 1 describes the electronic structure 
of the a = I (left) and r (right) electrodes and the coupling between them and the ladder. The 
specific electronic properties of the electrodes will not appear explicitly in our treatment; they are 
considered only as sinks of charges and will be included in the calculations via phenomenologi- 
cal parameters. The electrode-ladder interaction, given by the tka.sj matrix elements is in general 
not limited to the ladder sites nearest to the electrodes but it can have a longer range. This will 
qualitatively account for the fact that the real structural molecular arrangement in the active region 
does not in general have a linear chain topology —it is expected to rather be a complex network- 
so that interactions do not necessarily have to be confined to the electrodes' nearest neighboring 
molecular blocks. 

The Ansatz for £j(V) allows to account qualitatively for two different physical scenarios: with 
increasing V > 0, the built-in field K can be progressively counter-acted by the voltage-dependent 
term (j/N)eV, the energy levels become increasingly aligned and the current can dramatically 
grow due to the reduction of tunnel barriers between the sites. For negative bias, transport is 
blocked up to a critical bias —V*; the energy scale \eV* | is roughly related to the situation where the 
//-level at a site j and the L- level at site j + 1 become nearly resonant (additionally, at least the first 
site of the //-strand should move above the left chemical potential in order to have a current flow). 
Once this happens, an exponential current onset can occur. We have used two different modeling 
strategies based on two possible different microscopic transport mechanisms: hopping transport 
(HT) and coherent tunneling (CT). HT seems to be reasonable, since it is expected that the active 
device region displays some degree of static and dynamic disorder which may suppress coherent 
tunneling pathways over long distances. We have used a reduced density matrix approach 22 (for 
HT) and the Landauer formalism 23 to compute the I-V characteristics of the model in the CT 
regime (see the Supporting Information for details of the calculation). 

The results (scaled by an arbitrary factor) are displayed in [figure] [4] [] 4 in comparison with the 
experiments. It is obvious that the weak (strong) thickness dependence for positive (negative) volt- 
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Figure 4: Calculated I-V curves in the (a) hopping and (b) coherent transport regimes compared 
with the measured characteristics. The inset shows a schematic drawing of the two legs electronic 
ladder used to model charge transport in the device. The basic model parameters are: (a) ?o=150 
meV, t^ L =50 meV, t% ~ L =0, rf = rf = rf = rf=25 meV, and K=0.5, 0.37 and 0.27 eV for N=3, 
4, and 5; (b) f =150 meV, rf~ L =150 meV, fJ^ L =50 meV, rf = rf = rf = rf=100 meV, j8=1.2, 
and K=0.66, 0.6 and 0.55 eV for N=6, 7, and 8. For a definition of the model parameters we also 
refer the reader to the Supporting Information. 
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ages can be qualitatively described within both transport regimes. Notice, however, the different 
number of sites needed in both regimes to account for the experimental curves, This is related to the 
fact that tunneling processes can still be efficient for short lengths, so that almost no zero-current 
gap is obtained. Only with increasing length tunneling transport is blocked at low bias. 

Some basic features of the transport process can be summarized as follows: for positive volt- 
ages the onset of current takes place approximately at the critical bias where the first site of the 
HOMO-strand moves above the left chemical potential /!/ (which is kept fixed at \i\ = Ep while 
\x r = Ep + eV). Under this condition, charge transport can take place. A further increase of the bias 
leads to an increase of the current till all the electronic levels along the ladder become aligned (the 
external potential cancels the built-in field k). This provides the maximum current in this model; 
for larger positive voltages the misalignment of the levels develops again and negative differential 
resistance can eventually be obtained. For negative voltages, the current is almost zero in the case 
of coherent tunneling, since this process is very sensitive to the presence of energy gaps along the 
ladder. This behavior can be lifted by allowing for long-range interactions between the electrodes 
and the ladder states (see Eq. 7 in the Supporting Information). On the contrary, hopping transport 
can also take place in the negative bias regime with only nearest- neighbor interactions. 

Our results suggest that the effective microscopic transport mechanism may include both HT 
and CT contributions, each of them becoming dominant in different voltage regions. This may be 
the reason for the different slopes in the /(V = const. ,T) vs. T 1 plots for positive and negative 
voltages ([figure] [2] []2). A pure CT model cannot account for a T-dependence; HT contributions 
can however explain it. In the HT regime, a different qualitative behavior for positive and negative 
voltages is expected on the basis of a simple argument: the nearest neighbor inter-site hopping rates 
within the ladder include (see Eq. 4 in the Supporting Information) thermal factors exp(— (ej +1 — 
£j)/ksT), whenever ej +1 — e p - > (s,p = H,L). For s ^ p (sites on different strands), one obtains 
£j + l — £j=AH-L—K + eV/N; from here it is obvious that for V > the temperature dependence is 
mainly determined by the gap A H _ L <C kgT (large slope). For negative bias, A H _ L can be partially 
compensated by the term — K — e\V\/N and thus the temperature dependence may be weakened. 
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A more accurate quantitative analysis requires a systematic experimental study of the temperature 
dependence, taking also the influence of the transport layers into account. 

Conclusions. In this Letter, we unambiguously observe signatures of Zener tunneling in a 
organic molecular semiconductor: The reverse currents in organic pm-diodes show geometry, field 
and temperature dependence in excellent agreement with this tunneling effect not observed so far in 
such systems. The realized pm-Zener diodes have excellent device properties and allow important 
application possibilities, e.g. to reduce crosstalk in passive matrix arrays. The basic effect of 
valence- to conduction-level tunneling could be used in further device concepts, such as devices 
employing tunable resonant tunneling. Using a minimal model Hamiltonian, we can describe 
semi-quantitatively the measured electrical response. However, both limiting cases investigated 
here, incoherent hopping or coherent tunneling, can describe the observed experimental features. It 
therefore seems appealing to expect that the dominant charge transport mechanism may sensitively 
depend on the sign of the applied bias. A further point to be considered is the coupling to dynamical 
degrees of freedom which may induce polaron formation, so that the transport process will be 
rather mediated by dressed electrons or holes. We have not taken these issues in our model which 
is only addressing electronic degrees of freedom. Further work in this direction including detailed 
first-principle electronic structure calculations is thus needed. 
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Mott-Schottky-curves for the free available materials. Furthermore we present a more detailed 
description of the microscopic model applied to describe the charge transport characteristics in the 
organic Zener diode. 

This material is available free of charge via the Internet at http : / /pubs . acs . org/. 
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Experimental Methods 

Preparation. The devices are prepared on cleaned (supersonic treatment in acetone, ethanol and 
isopropanol for 5 minutes and 02 plasma etching) glass substrate by thermal evaporation of metals 
and organic materials under high- vacuum conditions (base pressure < 10~ 5 Pa) without breaking 
the vacuum. The layer evaporation is done in the following sequence: Aluminium (bottom con- 
tact, thickness 50 nm), HTL (50 nm), IIL, ETL (50 nm), Aluminium (top contact 50 nm). The 

thicknesses are confirmed by quartz crystal monitoring and the active area is 6.38 mm 2 . 
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The main materials used are: MeO-TPD: N,N,N9,N9 tetrakis (4-methoxyphenyl)- benzidine. 
NPB: N,N9- di(naphthalen-l-yl)- N,N9)-diphenyl- benzidine. BAlq2: aluminium- (Ill)bis (2- 
methyl-8-quninolinato)- phenylphenolate. Bphen: 4,7-diphenyl-l,10-phenanthroline. RE68: Tris 
(1-phenylisoquinoline) iridium(III). For the p-type doping, of the HTL the material NDP2 (No- 
vated AG, Dresden) can be used as well as the material F4TCNQ: 2,3,5,6- tetrafluoro-7,7,8,8- 
tetracyanoquino- dimethane as freely available p-dopant giving equivalent electrical properties (see 
[figure] [1] [] 1). The n-type doping of the ETL can be achieved by co-evaporation of Bphen and pure 
Caesium. Likewise, the material NDN1 (Novated AG, Dresden) acts as an effective n-type dopant 
(see [figure][l][]l). 

For the electronic characterization a Keithley 2400 SMU is used for I-V investigations while a 
HP 4284 LO?-meter is used for the C-V and C-f studies. The I-V curves are done with a voltage 
sweep (from -10V to 4V and back) and the C-V measurements are prepared at 1 kHz with an 
alternating sine signal (RMS value) of 20 mV. 

Organic Zener diode with different molecular dopants. Owing to the highly efficient doping of 
the hole and electron transport layers, the organic Zener diode becomes a reproducible and stable 
device with a low series resistance. Also, doping provides a second independent way to control the 
reverse behaviour. We can show (see Figure 1), that the effect of the reverse breakdown appears in- 
dependently of the used dopant. However, due to the different efficiencies of the dopant materials, 
one can observe an influence of the used dopants on the reverse direction. The number of charged 
dopant states at the interface between the doped layers and the intrinsic interlayer determines the 
thickness of the depletion region which controls the reverse current. Accordingly, one can further 
extract from Figure 1 that the molecular doping of MeO-TPD by F4TCNQ is less efficient than by 
NDP2 (for the same mass-ratio,corresponding to a molar percentage of approximately 8mol% for 
F4TCNQ and 6mol% for NDP2, respectively). Nevertheless, the series resistance of the transport 
layers is not significantly affected and hence no change in the forward direction is obtained. 

C-f and C-V investigations. In the frequency range from 50Hz to 1MHz, the organic Zener 
diode provides a phase and impedance response (see Figure 2 a) which fits to a simple RC-unit 
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Figure 1: I-V curves for different p-type and n-type dopants. The device consists of an intrinsic 
interlayer (7nm) of Balq2:NPB sandwiched between a hole and an electron transport layer (both 
50nm) of MeO-TPD and Bphen, respectively. The inset shows the I-V curve of a commercial 
inorganic silicon Zener diode with a breakdown voltage of 2.7V. 

as the equivalent circuit. Since the doped transport layers have a small series resistance and also 
a small geometrical capacitance (thickness 50nm) the RC-unit belongs to the intrinsic interlayer, 
which has the dominant RC-time in this stack (inset of Figure 2 a). One can thus identify the 
capacitance with the geometrical thickness of the depleted area. The series resistance describes 
Ohmic losses along the contacts and the transport layers, while the resistance parallel to the ca- 
pacitor models the influence of charge recombination and leakage currents for the reverse biased 
junction. This correlation between the capacitance and the thickness of the depletion, which is 
the sum of the interlayer thickness and the interface depletion, can be seen in Figure 2 b. For a 
reverse voltage between OV and -6V we obtain a linear dependence in the Mott- Schottky plot. 
This means, the C-V curves follow 



1/C 2 



2 Na + Nd 
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(V-V bi )+df. 



(1) 



Here, and are the densities of charged dopant states, d ; is the IIL thickness, A the active 
device area, V the applied voltage, and V/„- the built-in voltage of the diode. Accordingly, there is 
a voltage dependent depletion at the interfaces between the transport layers and the intrinsic inter- 




0.1 1 10 100 1000 
Frequency [kHz] 

(b) 

0.008r , , , , , , , ,— 




Voltage [V] 

Figure 2: C-f and C-V curves of the organic Zener diode, (a) C-f measurements of the or- 
ganic Zener diode. In accordance to the frequency dependence of the phase and the impedance 
one can assume a simple RC-unit as the equivalent circuit (inset), (b) Mott-Schottky plot of 
the device shown in (a). The device consists of Al(lOOnm)/ MeO-TPD:F4TCNQ[4wt%](50nm)/ 
Balq2:NPB[50wt%](7nm)/ Bphen:Cs(50nm)/ Al(100 nm) 
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layer as predicted by the Mott-Schottky rule. This means that the assumption of two rectangular 
space charge regions separated by an intrinsic layer is fulfilled in our pin-structure. With respect 
to this, we can expect a linear voltage drop over the intrinsic layer. It should be further mentioned 
that the Mott-Schottky rule is not valid for the forward biased device and also in the range of high 
negative voltages. The contributions of drift, diffusion current and recombination processes, which 
are not taken into account in this simple model, cause a collapse of the assumed equivalent circuit. 

Theoretical Methods 

Using the model introduced in the main text, Eq. 1, we can formally treat different transport scenar- 
ios. In the following we will address two limiting cases: incoherent, thermally assisted transport 
and coherent tunneling. 

Incoherent transport. The model in Eq. 1 has been applied to calculate the I-V characteristics 
for different lengths N. We may expect that static and dynamical disorder in the active region 
will make a fully coherent transport scenario very unlikely. Only for narrow regions (short length 
scales) quasi-coherent transport might still be efficient, though this may also depend on the voltage 
region. The electrical response of the ladder will thus be first studied in the limit of complete 
coherence loss. For this, a reduced density matrix approach 1 has been used and only its diagonal 
components (site populations) will be considered. Off-diagonal terms are assumed to have decayed 
on any relevant time scale. As a result, a set of classical rate equations can be derived (setting in 
this case the next nearest neighbor coupling tg ~ L to zero, but keeping t^~ L ^ 0): 

jPnAt) = T hn + T 2,n (f2{^)-PH,n{t)) + rHH,n+lPH,n+l{t) + rHH,n-lPH,n-l{t) 

+ 7hL,n-lPt,n-l(t) + JhL,n+lPt,n+u(t) ~ (YhHm-1 + YtIH,n+l + YtiL,n- \ + YHL,n+l) PfJ ,n(t) 

J t P L,n(<) = n,„ (/l(e M L ) -Pl,„(0) +r 2 ,„ (/ 2 (e M L ) -P L ,n{t)) + YLL,n+lPL,„+l(t) + YLL,n-lPL,n-l(t) 

+ YLH,n-lPH,n-l(t) + YLH,n+lPH,n+ll(t) ~ {YLL,n-\ + YhL,n+\ + YLH ,n-l + YLH ,n+l) PL,n{t) ■ 

(2) 
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Here, PH,n(t) and PL,n{t) are the corresponding populations and the transition rates y can be 
expressed in terms of the tight-binding hopping parameters along the ladder of the model Hamilto- 
nian in Eq. 1 (main text). f s= i^{e n ' ) are equilibrium Fermi functions of the electrodes evaluated 
at energies e n ' . The parameters T sn = T s ^ n=m describe phenomenologically the charge hopping 
rates into and out of the electrodes and can be formally expressed via the coupling terms tka,sj us " 
ing Fermi's Golden Rule. Assuming them to be energy-independent constants within the so called 
wide-band limit, 2 we make the following Ansatz: 

^l,nm = SniSmitf, ^ r,nm = S n=N d m=N T^ if n G H strand 

r/,wn = S n ^ N+ i8 m=N+ iTf, T rinm = 8 n= 2N8 m= 2N^r if e L strand. (3) 

For the hopping rates we further assume: 3 

e S p(V)-<&(V) 

Yspj = {fPfe- w if e s p (V)-eP(V)>0, 

Yspj = (n 2 if e s p (v)-e?(V)<0, (4) 

where s and p denote the H or L strands, while n,m are site indexes along the strands or between 
them. Depending on the energy difference between the sites there may be an energetic barrier to 
climb or an activation-less hopping process. In the stationary regime, dPn.nl dt = dPL,n/dt = 0, 
and a linear set of equations for the populations can be obtained. Finally, using the definition of 
the stationary electrical current coming from lead a, I a (V) = — e < dN a /dt >, one finally gets 

ho,(V) = e -(r^[f r (elv)-p H , N (v)]+r^[f r (e^v)-p L , N (v)}) 

- e -(rf [f l {etv)-PiAy)]+r\[f l {e[,v)-p^ L {v)])- (5) 

Coherent tunneling. To explore to which degree coherent charge propagation could contribute to 
the observed length dependence of the I-V characteristics, we have also used the Landauer theory 2 
for the model of Eq. 1 . Hereby, the key quantity to be computed is the energy dependent quantum 
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mechanical transmission probability T{E), which encodes the electronic structure of the system as 
well as the coupling to the electrodes. T{E) can be computed according to the standard expression 
T{E) = Tr[G a (E)T r G r (E)Ti]. G r,a {E) denote retarded (advanced) ladder Green functions includ- 
ing the interaction with the electrodes. The electrode-system coupling matrices T^ r are defined 
below in 7. The Green functions can be calculated via Dyson's equation. 4 The electrical current 
through the system is then obtained (with /!/ = Ep and \l r = Ep + eV being the chemical potentials 
of the electrodes) as: 

I<y) = jfdE\f(E-in)-f(E-Hr)]T(E). (6) 

However, in order to have a non-zero current for negative voltages, it turned out that coupling 
beyond nearest neighbors, i.e. t#~ L ^ 0, as well as a long-range coupling T^ r to the electrodes was 
required. The latter was now modeled as: 

r /;WM = 8 nm r? eP n , r,„ m = 5„ w rf ifneH strand 

r /;WW = 5n m rfe^ n , r r ,„ m = 8 nm T L r e^ N -^ if n e L strand. (7) 

The reason for the strong current suppression at negative biases is related to the fact that in a 
Stark-ladder (linear slope in the onsite energy profile along the chain) as that used in the present 
calculations all electronic states are localized in the thermodynamic limit, so that for negative 
bias (where the slope of the ladder is very large) charge tunneling is essentially blocked. The 
calculated I-V curves are shown in the main text, see Fig. 4. The theoretical curves were scaled 
by an arbitrary factor, since the experiments measure current densities and in our approach only 
currents are calculated. Hence, a direct quantitative comparison between theory and experiments 
turns out to be very difficult in general. Nevertheless, the trends in the length dependence can also 
be approximately described by using this coherent transport model. 

In general, it may also be possible that both coherent and incoherent pathways are competing 
on different bias regimes. To clear this point requires however a more detailed microscopic study 
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which should include a first-principle estimation of the different parameters used in the model as 
well as the inclusion of a realistic arrangement of the molecular units in the device region. 
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